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G^ ' Abstract 



PLi| a new method is described for determining the quantum state of correlated multimode radiation by 

interfering the modes and measuring the statistics of the superimposed fields in four-port balanced 
homodyne detection. The full information on the A^-mode quantum state is obtained by controlling 
^ I both the relative amplitudes and the phases of the modes, which simplifies the reconstruction of 

density matrices to only A'^ + 1 Fourier transforms. In particular, this method yields time-correlated 
multimode density matrices of optical pulses by superimposing the signal by a sequence of short 
local-oscillator pulses. 



Recently, increasing interest has been devoted to the problem of reconstruction of the quantum 
state of optical fields from measurable data. The feasibility of reconstruction of single-mode density 
matrices was first demonstrated by Smithey et al. using optical homodyne tomography. In the 
experiments they determined both the Wigner function and the density matrix in a field-strength 
basis. Later, the problem of direct reconstruction of the density matrix in a field-strength basis and 
the photon-number basis has been studied in Refs. and respectively. All these methods are 
based on the fact that in balanced homodyne four-port detection the field-strength distributions 
of a signal mode can be measured and knowledge of the field-strength distributions for all phases 
within a vr-interval is equivalent to knowledge of the quantum state of the mode . Whereas for 
determining the Wigner function a three-fold integral must be calculated, the direct reconstruction 
of the density matrix elements can be accomplished with two integrals. The method has first been 
extended to two-mode quantum states ||5| and later on the general situation for reconstructing 
A^-mode quantum states from A^-fold joint difference-count distributions has been investigated 
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1^]. In these schemes the determinations of A^-mode density matrices and phase space functions, 
respectively, require 2N- and 3A^-fold integral transforms of the measured data. 

In the present paper we propose a new method for determining A^-mode density matrices [N 
> 1) from the difference count distributions measurable in balanced homodyne four-port detection 
with pre-superimposed signal fields. Apart from the advantage that only two photodetectors are 
needed in that scheme, the reconstruction of the density matrices only requires A^+l rather than 2N 
integral transforms. The method does not only enable one to measure the (joint) density matrices 
of correlated A^-mode fields very efficiently, but also yields A^-mode density matrices of optical 
pulses in terms of time-localized nonmonochromatic modes. The basic quantities in our scheme are 
the distributions of the sum field strengths of two or more modes undergoing both amplitude and 
phase control. It is shown that these distributions yield the full information on the quantum state, 
provided that they are recorded as functions of the relative field amplitudes of the modes and for 
phases within 7r-intervals. 

We will illustrate the principle by dealing with the two-mode case. An extension of the corre- 
sponding results to A^ modes (A^ > 2) is straightforward. Consider two typical experimental schemes 
as given in Fig. 1. First, let us assume that two (correlated) modes can be used separately as input 




Figure 1: Two possible schemes for reconstructing two- mode density matrices, 
(a) Two modes (ai and 02) are mixed by a beam splitter and one of the interfering output modes 
(6) is used as signal mode in balanced homodyning in order to measure the sum field strengths of 
the two modes (LO, strong local oscillator). 

(5) A signal pulse and a sequence of two short (strong) local-oscillator pulses with envelopes fi{t) 
and /2(t) are superimposed in balanced homodyne detection. 

fields in an interferometer, such as a beam splitter. The superimposed light in one output channel 
of the interferometer is used as a signal mode in balanced four-port homodyne detection [Fig. 0(a)]. 
Second, two-mode density matrices of optical pulses can be measured [Fig. |^(6)] by superimposing 
a signal pulse with a sequence of two short (strong) local oscillator (LO) pulses and recording 
the (time-integrated) difference-count distributions Q. The time-localized (nonmonochromatic) 
modes are defined by the short local-oscillator pulse envelopes fi{t) and f2{t), and the measured 
difference-count distributions yield the distributions of the sum field strengths of the two modes 



[10|. In both schemes, the relative amplitude of the two modes and their phases are controlled. 
Consider a two-mode optical field with photon destruction operators (A; = 1, 2) and introduce 
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the scaled field-strength operators of the modes as 

Ffc(^fc) = |F|(afce-*'^'=+4e^^'^), (1) 

and Fk{fk)\^k,^k) = ^k\^kifk)- The mode amplitude |F| specifies the kinds of fields (electric, 
magnetic) under study, to compare with experiments it can be related to the shot noise [0]. In the 
field strength basis the density matrix is given by 

= (^) j dy^ j dy2e-'^y^^^+y^^-'>^{zi,Z2,^i,^2), (2) 

where 

Zk = Zk{yk,J"k) = ^yk' + J"^/\F\\ (3) 

V'fc = ipk{yk,J^'k) = ^k- arccot[yfc|Fp/J^fc] (4) 

{J-'f. > 0) . The characteristic function ^' in Eq. (^) is closely related to the average of the two- mode 
coherent displacement operator and can be written as 

'^{zi,Z2,ipi,'>p2) = (exp iziFi{ipi) +222^2(^2) )• (5) 

We see that ^' as a function of the Zk (for given ipk) is nothing but the characteristic function of 
the joint field-strength distribution 

= (^)'/ / d^2e-^(^i^i+^2^2)^(zi, 22,^1,^/^2). (6) 

Thus, when the joint probability distributions of the two field strengths are measured for all values 
of the phases ipk within vr intervals, two- fold Fourier transforms yield their characteristic functions 
and two more Fourier integrals must be performed to obtain the density matrix in Eq. (^). 

Let us now suppose that in place of the above mentioned joint field-strength distributions the 
distributions Ps{J^, a, V'l, ^^2) of the sum field strength 

F = Fi(V'i) cosa -|- -F2(V'2) sina (7) 

are measured, where aG(0, ^vr) controls the relative field amplitude. The probability distributions 
Ps{J^,a,'tpi,ip2) and pj(^i, .F2, V'l, V'2) are related to each other as 

Ps{T,a,il'i,ip2) = J dJ^i J dJ^2PiiJ^i,^2,i'i,tp2)SiJ^-J^icosa-J^2sma), (8) 

which together with Eq. (^ implies that 

Ps(^, a, V'l, "02) = — /d2;e~*^'^^'(2;cosa,2:sinQ;,'i/;i,'!/'2), (9) 
zvr J 



and hence 



where 



^(zi,Z2,Vi,^2) = yd.Fe'"^Ps(.^,a,Vi,V'2), (10) 



Z = Z{ZI,Z2) = Jzf + Z^, (11) 



3 



a = a{zi,Z2) = arctan(z2/-zi)- (12) 

This means that from measurements of the superimposed hght (^) for different a we can "to- 

mographically" reconstruct the joint characteristic function (0). Combining Eqs. (|2|) and (p!o|) 
yields 

= (^) / dyi j Ay2e-'^y^^^+y^^^^ J d^e^^-^p^l^, /?, V-i, V'2), (13) 
where [according to Eqs. (|3|), (11), and ([T^)] 



y ^ y(yi,y2,^{,^^) = Jyj + vl + (-^f + ^?)I\F\^ , (14) 



/? = /3(yi,y2,-^{,-?^2) = arctan W ^ 



(15) 



and ^/'fc is given in Eq. (^). Equation (13) exphcitly shows the feasibihty of reconstructing the 
two- mode density matrix from the measured sum field-strength distributions [V'A: ^i}Pk ~'^^^k)^ P 
G(0,i7r)]. 

Let us briefly comment on the effect of nonperfect detection. When the photodetectors have effi- 
ciencies 7] that are less than unity, the characteristic function of the (sum) field strength distribution 
in Eq. (|l3|) is replaced according to |^] 



where ps(-^5 /3i "01) "02! f?) denotes the sum field distribution measured with a detector of quantum 
efficiency rj. Consequently, Eq. (|l^) is generalized as 

1 \2 



— j I dyi I dy2e''^y'^'+y^^^hy 1^1 jdTe'y^ps{J',P,^i,^2;v)- 



(17) 



Processing data from a real experiment requires a careful handling of the noise to avoid troubles 
arising from the exponential of y"^ . 

Let us return to the measurement schemes proposed above. Equation (^^ yields the general 
relation of the measured distributions to the corresponding two-mode density matrices. In the 
first scheme [Fig. |l](a)] the difference count distributions directly yield the desired distributions 
p[T, (3, ipi,ip2',v) of the sum field of the two signal modes under study. The parameters ipi, 'ip2, and 
f3 are related to the input-output relations of the interferometer. Measurement of the difference- 
count distributions on a sufficiently dense grid of points {tpi,tp2, f3), with ^pk & ifk — T^^^k) and (5 G 
(0, ^vr), is then equivalent to measurement of the quantum state of the two- mode field. In the second 
scheme [Fig. |l](&)], two-mode density matrices of optical pulses can be measured by superimposing 
the signal pulse with two short (strong) local oscillator (LO) pulses of envelopes fi{t) and f2{t)- The 
desired information on the two-mode quantum state is obtained by recording the (time-integrated) 
difference-count distributions and performing the amplitude and phase control by controlling the 
relative strength and the positions of the LO pulses, for more details see [p. Note that this method 
is not only of interest for studying pulses, it also yields insight in the time-dependent correlation 
properties of stationary fields at the level of the full quantum statistical information. 
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In certain cases of pulse measurements, when the signal pulse and the LO pulses are produced 
by different sources, it may become impossible to achieve complete phase control. However, the 
difference of the LO phases, /\ip = ip2 — '(pi, can be controlled and we measure the phase averaged 
probability distribution, 

Ps{T,(3,Aij) = (27r)-i J d^ip,{T, P, Vi, V'l + A^-). (18) 

Concequently, we can reconstruct the phase-averaged (with respect to V'l) version of the density 
matrix in Eq. (0). Although losing some information, these averaged density matrices still contain 
interesting informations on the temporal correlations of the signal field. 

It is straightforward to extend the method to measurements of A^-mode density matrices (A^ > 
2). Again, we can measure the weighted sums of the field strengths in a homodyne setup with two 
detectors. The pre-superimposed signal can again be obtained using appropriate interferometric 
methods. In particular, in the scheme in Fig. ^(6) a train of A^ (strong) LO pulses must be 
used. The sum field strengths are now parameterized by A" — 1 relative intensities of the A^ modes 
under consideration, and the characteristic functions of their probability densities are related to 
the characteristic functions of the joint probability densities by generalizing Eq. (^). Then the 
reconstruction of the A^-mode density matrix can be accomplished with A^ + 1 Fourier integrals 
- one integral for obtaining the characteristic function of the density matrix from the measured 
data plus one integral per mode in order to determine the density matrix from its characteristic 
function. 

In conclusion, we have shown that the reconstruction of multimode density matrices of light 
can be performed by measuring amplitude and phase controlled distributions of the sum fields. 
This simplifies the reconstruction of multimode density matrices by reducing both the number 
of photodetectors needed and the number of the desired integral transforms. Applications of the 
method for determining the quantum states of correlated field modes and of time-correlated mul- 
timode quantum states are discussed. The latter allows to perform time-dependent correlation 
measurements yielding the full quantum statistical information. 

We thank R. Sauerbrey for useful discussions. This work was supported by the Deutsche 
Forschungsgemeinschaft . 
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